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Secondary Heavy Quark Production in Jets through Mass Modes 
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We present an effective field theory method to determine secondary massive quark effects in jet 
production taking the thrust distribution for e + e~ collisions in the dijet limit as a concrete example. 
The method is based on the field theoretic treatment of collinear and soft mass modes which have to 
be separated coherently from the collinear and ultrasoft modes related to massless quarks and gluons. 
For thrust the structure of the conceptual setup is closely related to the production of massive gauge 
bosons and involves four different effective field theories to describe all possible kinematic situations. 
The effective field theories are patched together and allow for a continuous description of infinitely 
heavy to arbitrarily small masses keeping the exact mass dependence of the singular terms. The 
mass mode field theory method we present here is in the spirit of the variable fermion number 
scheme originally proposed by Aivazis, Collins, Olness and Tung and can also be applied in hadron 
collisions. 



I. INTRODUCTION 

By now jet physics has reached a high level of preci- 
sion allowing for an accurate description of the strong 
interaction. Fundamental parameters of QCD such as 
the strong coupling constant as well as nonperturbativc 
properties of hadrons like parton distribution functions 
can be determined with a continuously improved accu- 
racy from high precision data samples. On the theo- 
retical side this has been possible thanks to high order 
loop calculations and the summation of large logarith- 
mic terms. Computations of jet cross sections for mass- 
less quarks belong to the well known and unambiguously 
defined exercises in perturbative QCD based on a num- 
ber of rigorous factorization proofs. On the other hand, 
as far as massive quarks are concerned it is fair to say 
that their treatment is not coherent throughout the lit- 
erature. Different schemes for massive quarks exist which 
differ in the resummation of logarithms and in the inclu- 
sion of formally subleading contributions. An approach 
capable of describing quark mass effects starting in prin- 
ciple from very small masses when the quarks are inside 
hadrons and stretching up to ultra-heavy masses in the 
decoupling limit was provided by Aivazis, Collins, Ol- 
ness and Tung (ACOT) [l|, 0. Their variable fermion 
number scheme is based on the separation of close-to- 
mass-shell and off-shell modes and allows to factorize 
infrared-safe hard coefficient corrections from low-energy 
parton distribution terms involving logarithmic mass ef- 
fects. In this respect the concept behind the ACOT 
scheme is along the lines of effective field theory methods 
such as the soft-collinear effective theory (SCET) frame- 
work and can be readily incorporated into it. As we 
show in this paper, the resulting effective theory frame- 
work can be based on the inclusion of collinear and soft 
"mass modes" 1 together with the existing collinear and 





1 This differs from the terminology used in Refs. [f| and [(J where 
only the soft massive modes were referred to as the mass modes. 



FIG. 1: Diagrams at O(oi) for virtual and real secondary 
radiation of massive quark pairs in primary massless quark 
production. 



soft massless partonic modes. While the collinear and 
soft massless partonic modes typically have different in- 
variant masses depending on the observable under con- 
sideration, the collinear and soft mass modes have the 
same typical invariant mass. This can lead to compli- 
cated patterns of scale hierarchies that might even vary 
substantially within a single distribution. 

In this and a subsequent paper we apply the mass mode 
SCET concept to describe the secondary production of 
massive quarks in the e + e~ thrust distribution. While we 
aim for the description of thrust, the described method 
is general and can be applied with possible adaptions 
to other processes as well. Recently, a factorization ap- 
proach was derived within SCET that can be applied for 
thrust induced by primary heavy quarks, i.e. by heavy 
quarks that are produced directly from the hard jet cur- 
rent [H, @ . For the description of primary heavy quarks 
the mass modes are involved as well, but they play a more 
simplistic role since massless collinear modes are directly 
tied to the massive collinear modes and soft mass modes 
turn out to only lead to virtual effects. So the result- 
ing situation is almost identical to the one in massless 
SCET. For secondary heavy quark production, i.e. for 
heavy quarks produced from gluon radiation off primary 
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massless quarks, Fig. [T] the collinear and soft massless 
as well as massive modes can all constitute independent 
degrees of freedom, each contributing to thresholds and 
dynamical effects. Here the dynamical mass effects enter 
in different ways depending on the relation of the mass 
scale to the kinematic scales relevant for massless quark 
production. Depending on the mass value the scale hi- 
erarchies can vary substantially requiring different types 
of effective field theory setups that have to be patched 
together in a connected way to allow for a continuous de- 
scription of the distribution. The secondary quark mass 
effects in thrust thus represent a non-trivial prototypical 
showcase for the mass mode concept. 

Concerning the numerical size the effects from sec- 
ondary heavy quarks in thrust are 0{a 2 s ) in fixed-order 
perturbation theory and certainly small if the hard scale 
Q is much bigger than the heavy quark mass. Never- 
theless the numerical effects involve logarithms related 
to renormalization evolution with variable flavor number 
already at the leading-logarithmic level and are impor- 
tant if the hard scale is only a few times larger than the 
mass. They are certainly essential for a precision anal- 
ysis of e + e~ event shapes for Q < 35 GeV concerning 
the description of bottom mass effects. So studying the 
secondary heavy quark effects in thrust is also of phe- 
nomenological interest. 

An interesting conceptual issue, which can also lead 
to considerable computational simplifications, is that 
the description of secondary massive quarks radiation is 
closely tied to radiation of massive gauge bosons. All 
non-trivial conceptual aspects and also some computa- 
tional issues can be conveniently addressed considering 
the setup with massive gauge bosons which we mostly 
concentrate on in this first paper. Details on the com- 
putation of the effects of secondary massive quarks and 
extensive numerical studies are postponed to Ref. [7|. 

This paper is organized as follows. In Sec. [TT] we ex- 
plain the outline and the scope of the mass mode method. 
We also show how secondary massive quark radiation is 
tied up to the description of massive gauge bosons. So 
the required effective field theory framework can be con- 
ventiently formulated concentrating on collinear and soft 
massive gauge boson radiation as a placeholder. Essen- 
tially all non-trivial conceptual and technical issues can 
be discussed and handled with massive gauge bosons. In 
Sec. IIIII we review the massless results for the thrust dif- 
ferential cross section. In Sec. IIV1 which is the central 
section of this paper, we explain our theoretical setup 
to describe corrections due to secondary massive parti- 
cles with mass modes and show the results for massive 
gauge bosons with vector coupling. Here we also ad- 
dress how the approach achieves a continuous description 
of the mass-dependence covering all possible hierarchies 
and we discuss why we do not miss any sublcading power 
corrections at the transition points between the different 
scenarios. Details on the actual calculations are given in 
Sec. |V]for the full theory results and in Sec. |VT]in the 
effective theory framework. As an outlook we describe 
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FIG. 2: Figure illustrating the dispersion method for the vac- 
uum polarization correction to the gluon propagator in the 
subtracted version with Il(q 2 = 0) = suitable for situations 
where the massive quark is not contributing to the renormal- 
ization group evolution. The explicit analytic form of the 
dispersion relations is discussed in Sec. IVHI 



briefly the two-loop computations required for secondary 
massive quarks in Sec. IVIII and show the impact of our 
results on the thrust distribution. Finally in Sec. lVIlIl wc 
summarize and conclude. 



II. THE OUTLINE 



We consider the effect of secondary massive quarks in 
the e + e _ thrust distribution for massless primary quark 
production, da hsht /dr. For thrust we use the definition 



i-r = i-£ 



|n -Pi\ 

E,l^-| 



, (1) 



where n is the thrust axis and the sum is performed over 
all final state particles with momenta pi and energies 
Ei? The most singular terms for small r (i.e. in the dijet 
limit) are governed by a factorization theorem that sep- 
arates the dynamical fluctuations at the center-of-mass 
energy Q (hard scale) , the typical invariant mass of each 
of the jets QX (jet scale) and the typical energy of large 
angle ultrasoft radiation QA 2 (ultrasoft scale), where the 
power counting parameter A ~ yfr in the tail region of 
the thrust distribution. For m ^> Q the effects of the 
massive quark decouple and the factorization theorem for 
the singular partonic cross-section adopts the well known 
form for massless quarks [5]. Schematically the form is 



dcr „ , „ 
dr 



(2) 



Here H denotes the hard contribution, J the jet func- 
tion describing collinear fluctuations, which is convoluted 
with the soft function S arising from soft large angle radi- 
ation. To avoid large logarithms each term is evaluated at 



2 Wc define the thrust variable r normalized with the c.m. energy 
Q, which is the sum of all energies and also agrees with the 
variable 2-jettiness §. For massless decay products this agrees 
with the common definition, which is normalized to the sum of 
momenta J2i \Pi\- 
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FIG. 3: Diagrams at 0{oi 2 a ) for virtual and real secondary 
radiation of gluons with mass M in primary massless quark 
production. 



its characteristic renormalization scale, so that pin ~ Q 
for the hard factor, fij ~ QX ~ <9v^" f° r the collinear 
jet function and (is ~ QA 2 ~ Qt for the soft function. 
Large logarithms between these scales are summed by 
the renormalization group factors which are implied. 

In order to determine the effects of secondary radia- 
tion of heavy quarks through gluon splitting (Fig. [TJ one 
has to deal with two issues: (i) the separation of dynam- 
ical modes depending on the relation of the mass m with 
respect to the hard, jet and ultrasoft scales and (ii) the 
calculation of the secondary massive quark contributions 
in connection with the appropriate number of flavors con- 
tributing to the renormalization group evolution. For jet 
observables and quantities depending only on the invari- 
ant mass of the secondary fermion pair the two issues 
can be completely disentangled. This is achieved by us- 
ing the fact that the quark pair polarization correction to 
the gluon propagator can be expressed as a dispersion in- 
tegral of a massive gluon propagator with the absorptive 
part of the quark- vector current correlator (see Fig. [2J . 
Using this dispersion integral method we can set 

up the conceptual formalism to describe the secondary 
massive quark pair radiation and to separate the mas- 
sive modes in connection to the massless collinear and 
soft modes by first considering real and virtual radiation 
of "gluons" with mass M, as shown in Fig. [3] All con- 
ceptual issues concerning the separation of modes, the 
determination of matching conditions and of matrix ele- 
ments within the context of SCET can be conveniently 
discussed with these massive gauge bosons. 

One has to consider collinear (for each jet direction) as 
well as soft massive "gluon" modes, all of which have the 
same mass and typical virtuality of order M. This dif- 
fers from the virtualities of collinear and ultrasoft mass- 
less modes which are of order QX and QX 2 , respectively, 
with QX ^> QX 2 . After having carried out all 0(a s ) 
calculations with the massive "gluons" the 0(a 2 ) effects 
from secondary massive quark radiation are obtained by 
the dispersion integration over the gauge boson mass ac- 
cording to Fig. [2] The dispersion integration can be car- 
ried out either in the unsubtracted or subtracted versions 
depending on whether the massive quark flavor gives a 



UV divergence and contributes to the RG-evolution or 
whether it is integrated out and does not contribute to 
the RG-evolution. 

In this paper we mostly concentrate on the conceptual 
aspects related to the collinear and soft gauge bosons 
with mass M and on their 0{a s ) effects in thrust. Wc 
aim at a continuous description for all mass values from 
M 3> Q, where decoupling takes place, down to the limit 
M — » 0, where the formalism continuously interpolates 
to the well-known case of a massless gluon. The linearity 
of the dispersion integration entails that the continuity 
is carried over also to the contributions coming from the 
secondary massive quarks. The continuity of the theo- 
retical description as a function of Q, the mass and the 
external kinematical variables (which is thrust r in our 
concrete application) is particularly important, since dif- 
ferent hierarchies between the mass scale and the jet and 
ultrasoft scales can arise within a single spectrum when 
the kinematic variables are varied within their allowed 
kinematic ranges. 

We stress that the mass mode method we present here 
is not tied to separating the problem into a treatment of 
massive gauge bosons and dispersion integrations. Once 
the field theoretic setup is identified (where the use of 
massive gauge bosons might certainly serve as a useful 
guideline) the approach can be formulated right away in 
terms of massive quark modes. In general this can be 
even necessary when jet observables are considered that 
do not depend in a universal way on the invariant mass 
of the quark pair. However, in cases where the dispersion 
method can be applied the results for the massive gauge 
bosons are universal and results for other types of sec- 
ondary massive particles (scalars, vector bosons) can be 
readily obtained simply by using the respective version of 
the dispersion relation. Moreover in many non-universal 
cases the computationally easier dispersion method can 
be used to at least correctly determine ultra-violet (UV) 
divergent and/or infrared (IR) singular terms, such that 
the remaining finite terms do not need to be carried out 
with regulators. For thrust this is the case concerning 
large angle soft radiation. 

An interesting technical aspect of computations with 
mass modes is that usual dimensional regularization is 
not capable of separating contributions from the collinear 
and soft sectors, see e.g. Refs. [I3 - ll4j |. This leads to 
unregularized divergences in the respective collinear and 
soft diagrams which appear to be related to the fact that 
dimensional regularization regularizes with respect to vir- 
tualities, whereas collinear and soft mass modes are sep- 
arated only by a boost. As we demonstrate in this work, 
these singularities are only spurious and drop out of prop- 
erly defined collinear and soft matrix elements. In par- 
ticular they are not related to any renormalization effect. 
Using the mass mode massive gauge bosons and the dis- 
persion method allows for a simple computational scheme 
to deal with these spurious singularities because the sin- 
gular terms can already be canceled at the level of the 
0(a s ) diagrams with the massive gauge bosons. This can 
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also lead to substantial technical simplifications in cases 
where the secondary massive quarks enter an observable 
in a non-universal way. 

We note that our effective theory setup involving the 
mass mode gauge bosons entails that they are distinct 
degrees of freedom and exist in parallel to the common 
massless collinear and ultrasoft gluons. While we primar- 
ily give them the role of placeholders for the mass mode 
collinear quarks (and the collinear gluonic modes they 
can interact with) the calculations we carry out in this 
work can be of direct interest by themselves and might 
well serve for concrete applications for example related 
to the electroweak theory. To be specific the mass mode 
gauge bosons might be considered within the context of 
an additional spontaneously broken SU(2) gauge theory 
where all SU(2) gauge bosons get a common mass M 
from a Higgs in the fundamental representation. Similar 
to Refs. [H, [ID, UU we write the corresponding SU(2) 
group theory factors using Cp, Ca and Tp and denote 
the gauge coupling as g with a s = g 2 /47r. This agrees 
with our notation for the QCD results and facilitates the 
interpretation of the results for the dispersion integra- 
tion used to implement the corrections due to the gluon 
splitting. 

As a final introductory remark we note that for the 
calculations of the secondary massive quark effects all 
computations are performed in dimensional regulariza- 
tion. Since in this paper we mostly focus on the concep- 
tual field theory and computational issues related to the 
massive gauge bosons we present all results in the limit 
d = 4 — 2e — > 4 to facilitate the discussions. Details on 
the d dimensional results together with the corrections 
caused by the secondary massive quarks are treated in a 
subsequent publication Q. 



III. THE MASSLESS FACTORIZATION 
THEOREM 

In this section we briefly review the known massless 
thrust factorization theorem to set up our notations and 
collect the perturbative QCD results at one loop. The in- 
clusion of the mass modes modifies the form of the mass- 
less factorization theorem, and the various mass mode 
contributions to the factorization theorem interpolate 
continuously to the respective massless contributions in 
the limit M — > 0. The factorization theorem reads |5j 

^ =(T Q 2 H (Q,fiH)U^\Q,fj,H,fJ-) 

x ds Ids' Jo(s',)j,j)Uy(s — s',^,fjbj) 



J des (Qr--^-e,ns)uP 



(^,M,Mfif)) ( 3 ) 



where ctq denotes the total partonic e + e~ cross-section 
at tree-level, H is the hard current matching condition, 
Jo the thrust jet function and So the thrust soft function. 



Large logarithms between the characteristic scales of each 
sector, fin , /ij and fig, and the final renormalization scale 

\i are summed by the evolution factors U# , U^p and 

U s satisfying the renormalization group equations, 



(4) 



d f 

V-^U ( j\s,(j,,hj) = J ds'jj(s - s' ',n)Uy\a' , 

(5) 



dfi 



(6) 



The superscript (1) in the evolution factors as well as the 
subscript in the matrix elements indicate an expres- 
sion in the massless theory without mass modes. The 
choice of fi is arbitrary, and the dependence on /i cancels 
exactly working at a particular order. In this work we 
adopt the choice /i = /is, such that the evolution factor 
Ug {£, p,s, fJ>s) — S(£) and can be dropped from Eq. 
The fact that other choices for /i can be implemented 
leads to interesting consistency conditions between the 
renormalization group factors which have been discussed 
in detail in Ref. [5(. We stress that beginning at 0(a 2 s ) 
all ingredients of the factorization theorem depend on the 
number of massless quarks n/. Currently Hq, Jq and the 
partonic contributions to So are known up to 0(a^_) [lQ- 
|20 ] and the anomalous dimensions up to 0{a 3 s ) [la. l21j. 

The hard function is Hq(Q,h) = \Cq(Q,[i)\ 2 with 
the massless Wilson coefficient Co(Q,/i) from matching 
SCET to QCD. At 0(a s ) it reads (with Q 2 = Q 2 + iO) 



C (Q,m) = 1 + 



ol s Cf 

47T 



-In 2 



-Q 2 



3 In 



7T 

- 8 + T 



(7) 



The current divergences are canceled by the current 
renormalization factor 



Z c (Q,fi) = l 



a s C F 

47T 



In 



z9l 

M 2 



(8) 



which yields for the leading order anomalous dimension 



1h{Q,h) = -Z c (Q,(i)fj,— Zc(Q,Li) + c.c. 



47r V V A* 



12 



(9) 



The thrust jet function Jq(s,[x) is obtained through a 
convolution of the hemisphere jet functions [6] for the 
two jet directions, 



J (s,fi)= / ds' Jo !n (s',fi)J ,n(s - s',fi) . 



(10) 
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The jet functions J n and J„ are matrix elements of 
collinear fields in SCET and describe the dynamics of 
the collinear degrees of freedom. The n— collinear jet 
function is defined as 
1 



4nN c Q 



xlm 



i / d 4 x e ir "- x (0\T{xn, Q mXn(x)}\0) 



(11) 



where the jet fields Xn{x),Xn,Q(0) denote quark fields 
multiplied by collinear Wilson lines and the invariant 
mass is r 2 ~ Qr+ . All color and spin indices are traced 
implicitly. For further details and the definition for 
4(Qr^,/i) we refer to Ref. @. 

The renormalized expression for Jo(s,/x) at one loop 
order reads 



H 2 J (s,n 2 ) = S(s) 



a.C 

47T 



' 1 5(s) (14 - 2tt 2 ) - 6 
s) In s 



0(s) 



+8 



(12) 



with s = s/fi 2 . The one- loop renormalization factor for 
the massless jet function reads 



6(s) 



H Zj(s,n) = S(s) + ^— ( 5{s) I + - I - - 

(13) 

which yields for the leading order anomalous dimension 



n 2 ^j(s,n) = - / ds' Z j 1 (s - s' , n) u— Zj(s' , fi) 



47T 



-16 



6(s) 




The thrust soft function Sq (£, n) describes ultrasoft radi- 
ation between the two jets and is based on a hemisphere 
definition, 

S (l, M ) = i J2 S ( £ ~ n ■ k * ~ n ■ k s) (0|F fi F„(0)|X s ) 



x (X s \Y^Yl(0)\0) 



(15) 



where kg (kg) is the momentum of the soft final state 
\X S ) in the right (left) hemisphere and Y n (x), Y n (x) are 
ultrasoft Wilson lines, i.e. 



Y n (x) = Pexp 
Y n (x) = Pexp 



ig dsn ■ A us (ns + x) 



-ig I dsn ■ A us (ns + x) 
Jo 



(16) 



In the tail region where the ultrasoft scale is larger than 
the hadronic scale Aqcd the thrust soft function factor- 
izes into a perturbative partonic part and a nonperturba- 
tive hadronic part. In the following we will consider just 



the partonic soft function. The renormalized expression 
at 0(a s ) is 



9{£)\n£ 



u S (i, n) = 5(1) + ^ j ^5(1) - 16 

(17) 

with £ = I j ' fi. Its renormalization factor at one loop reads 

[0(1)" 



,Zs(t) = S(i)-^ \± 2 5(i)-* e 



e 



, (18) 



+ 



and yields for the leading logarithmic anomalous dimen- 
sion 



M7s(AaO = - / d£'Zs X (t-£'^)^Z s (£',ti) 



d 

du 



a s C F 
Air 



16 



0(1) 



(19) 



The anomalous dimensions (to any loop order) satisfy the 
consistency relation 

1h(Q, h)5(s) + fJ?7j(s, n) + ^ls(s/Q, /i) = 0. (20) 



IV. THE MASS MODE SETUP 

When mass modes are included the mass scale can en- 
ter the field theory setup in addition to the hard, jet 
and ultrasoft scales. It is convenient to define the ratio 
Am = M/Q. There are different types of mass modes, 
collinear type modes and soft modes [lij. The n- and n- 
collinear mass modes have the scaling ~ Q(X 2 M , 1, Am) 
and ~ Q(1,A m ,Am)i respectively, while soft mass 
modes have ~ Q(Am,AmjAm)- The collinear mas- 
sive quark interactions are determined from a massive 
quark collinear Lagrangian 22j which is a straightforward 
generalization of the massless collinear Lagrangian. The 
soft mass mode gauge bosons couple to collinear quarks 
through soft Wilson lines [5j,|6| and the interactions of soft 
mass mode quarks and gluons among themselves is given 
by usual QCD interactions. An important aspect is that 
the mass mode collinear gauge bosons are defined with a 
soft mass mode-bin subtraction [l4[ in order to avoid dou- 
ble counting and maintain gauge invariance within each 
collinear sector. The collinear and the soft mass modes 
both have typical invariant masses \pp 2 n ~ \pP 2 s ~ 
and they must be considered together with the regular 
massless collinear and soft modes which have invariant 
masses of order QX and QX 2 , respectively, with A ~ y/r 
for thrust in the tail region. The collinear and the soft 
mass modes have the same typical invariant mass and are 
therefore just separated by boosts and rapidity. Depend- 
ing on the relative size of Q, Am and A, the mass modes 
have to be treated differently with respect to the mass- 
less collinear and ultrasoft modes. This has an impact 
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on the required matching and renormalization group cal- 
culations and therefore also on the form of the resulting 
factorization theorem. The different scenarios can be or- 
dered according to where the mass scale is situated with 
respect to the jet and ultrasoft scales, see Fig. [4] and [5] 
for graphical illustrations. Each scaling situation corre- 
sponds in principle to a different EFT setup. 

In the following we briefly summarize these different 
scenarios: 

I) M > Q > QX > QX 2 (A M > 1 > A > A 2 ): The 
mass M is larger than the hard scale Q and the 
mass modes are not contained in SCET, but inte- 
grated out when SCET is matched to the full the- 
ory. So the factorization theorem is the one for 
massless fermions up to the hard current match- 
ing coefficient which acquires an additional contri- 
bution due to the mass mode. The mass effects 
decouple for Am oo. 

II) Q > M > QX > QX 2 (1 > A M > A > A 2 ): In this 
scenario the mass M is in between the hard scale 
and the jet scale. The mass mode effects are vir- 
tual because the jet scale QA, the typical invariant 
mass for collinear particle radiation is below M. 
However, the mass modes contribute to the match- 
ing condition and the evolution of the current for 
scales above M. The mass modes are integrated out 
for scales below M. The factorization theorem is 
the one for the massless quarks concerning jet and 
soft function but receives mass corrections in the 
matching and evolution of the production current. 
The matching condition of the production current 
at the hard scale Q has the correct massless limit 
for A m -> 0. 

III) Q > QX > M > QX 2 (1 > A > Am > A 2 ): The 
mass is below the jet scale. Therefore massless 
and massive collinear modes belong to the same 
collinear sector both yielding real and virtual con- 
tributions to the jet function. The soft mass modes 
only arise through their virtual effects since their 
interactions with collinear modes induce invariant 
masses above the jet scale. They contribute to the 
current matching and evolution together with the 
virtual effects of the massive collinear modes ex- 
actly as in scenario II. Analogously to the previous 
scenario the massive modes are integrated out be- 
low the scale [Im ~ M leading to an additional 
collinear mass mode matching condition. The soft 
function agrees with the one of the massless factor- 
ization theorem. 

IV) Q > QX > QX 2 > M (1 > A > A 2 > \ M ): When 
the mass is below the ultrasoft scale the soft mass 
modes adopt the same invariant mass scaling as the 
massless soft degrees of freedom. However, they 
still have a mass and this must be taken into ac- 
count in the calculations. In this last scenario the 



massive modes are not integrated out and the fac- 
torization theorem has the form of the massless fac- 
torization theorem with the difference that all ma- 
trix element and evolution factors contain massless 
as well as mass mode contributions. 

An important aspect is that the results in each of 
these frameworks can be continuosly joined when vary- 
ing the mass or the relation between the hard, jet or soft 
scales. This is the case when scanning the thrust variable 
through the entire allowed spectrum. This yields a con- 
tinuous description valid for any scaling. At this point 
one might also worry whether the factorization theorems 
are subject to missing large power corrections at the point 
where two scenarios are patched together. We show be- 
low that this is not the case and that there is no generic 
loss of precision where the transition between different 
scenarios is carried out. 

We stress that the notation, the formulation of the fac- 
torization theorems and the organization of the RG evo- 
lution we employ is associated to the "top-down" evolu- 
tion, where the scale /z is equal to the ultrasoft scale fig. 
Thus only the current and the jet function evolution fac- 
tors Uh and Uj, respectively, appear in the factorization 
theorem. This also affects the interpretation and asso- 
ciation of the mass mode matching conditions discussed 
below. So in this RG setting the jet and the soft func- 
tions do not receive any massive quark effects when the 
jet and ultrasoft scales, respectively, are below the mass 
mode matching scale. Alternative ways to describe RG- 
evolution related to different choices for fi are possible 
and have been discussed in Ref. |6|. There are consis- 
tency conditions that relate the RG-evolution factors in 
the different schemes. Since the mass M appears as an 
additional scale in our mass mode setup, there are even 
more possibilities to describe RG-evolutions. Therefore 
the consistency conditions become more involved, as they 
also entail relations between mass mode matching condi- 
tions and renormalization group evolution with different 
anomalous dimensions and varying flavor number. As an 
example, if fi is set to the hard scale, the jet and soft 
functions have to be evolved upwards to larger scales. In 
this RG setup the jet and soft can pick up mass mode 
corrections in case the mass threshold is crossed, whereas 
the hard current does not receive any mass mode match- 
ing conditions. We discuss the impact of these scenarios 
in more detail in Ref. [7j. 

In the subsections IIV AllIV Dl we give a more detailed 
description of the formalism and the results for each of 
the scenarios for the application to thrust. We also ad- 
dress why the transition between the various scenarios 
is continuous and how this is related to the fact that no 
large power corrections are missed in the lower mass effec- 
tive theory description in the transition region. Overall 
this is conceptually based on the fact that the different 
effective theories are all matched to full QCD in an in- 
dependent way and should not be thought of as being a 
sequence of effective theories matched to each other. 
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FIG. 4: Localization of massless (ML) and mass modes (MM) in the p + -p~phase space according to their generic scaling for 
different hierarchies between A and Am- Modes with the same invariant mass are located on the same mass hyperbola. This is 
always the case for the collinear and soft mass modes. 



A. Scenario I: M > Q > QX > QX 



The massive gauge boson with vector coupling is inte- 
grated out at the hard scale Q, when SCET is matched 
to the full theory, yielding just a modification of the hard 
matching coefficient. It can be combined with the mass- 
less coefficient Co(Q, fi) to give a new Wilson coefficient, 
which reads at one loop order 



C I (Q,M,fx) = C (Q, f x) + SF m (M/Q). (21) 



The massive contribution SF m (M/Q) is computed from 
the full theory current form factor with a subtraction at 
(p + p') 2 = (i.e. in the on-shell scheme which ensures 
decoupling) since the mass modes do not contribute to 
the SCET renormalization group evolution. At one-loop 



order SF m (M/Q) reads [Tl|,[2 
a s C F 



5F m (x) 



, {(l + x 2 ) 2 [2U 2 (-x 2 )~ln 2 (-x 2 ) 

+21n(l + x 2 )ln(-x 2 ) - — 

3 

- (3 + 2a; 2 )ln(-a; 2 ) - 2x 2 - - 



(22) 



with x 2 = M 2 /(Q 2 + iO). In the limit M ->• oo the mass 
modes decouple, so SF m (x) — > for x — > oo. In the small 
mass limit x — ¥ (which is not supposed to be taken in 
this scenario), on the other hand, we obtain 



SF m (x) 



u s Cf 

47T 



ln 2 (-x 2 )+31n(-x 2 ) + ^ + 



2 

(23) 

which yields unresummed large logarithms. Thus in sce- 
nario I the correct massless limit cannot be obtained in 
the hard Wilson coefficient. Since the SCET setup in 
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II 



III 



IV 



M 



fi., ~ Q\ 



Ms ~ OA 2 



hard 



s e 
+ 



coll.ML 



MM 



hard 



coll.MM ■ 
soft MM 



= 3 

+ 



coll.ML 



soft ML 



soft ML 



hard 



+ 



coll ML 
icoll.M" 



: soft MM: 



soft ML 



hard 



coll.ML 
coll.M 



soft ML 
soft M 



fJ-M " 



MM = mass-mode, ML massless, M massive 



FIG. 5: The different scenarios depending on the hierarchy between the mass scale M and the hard, jet and ultrasoft scale. 
MM indicates mass mode scaling, ML the massless one. With M we denote modes that have a mass M but scale as their 
massless counterparts. The renormalization group evolution is also shown in the top-down evolution from the hard scale fin 
down to n = (is- When the mass scale is crossed the corresponding massive modes are integrated out. This leads to a matching 
condition and to a change in the evolution factor. In the case of secondary quark pairs production we evolve with nijght + iimass 
flavors above the mass scale and with n\i S ht below (n maS s is the number of quark flavors with mass M). 



this scenario is exactly the massless one, the factoriza- 
tion theorem adopts the form 

^ =Q 2 a Q \C I (Q,M^ H )\ 2 U i H ) {Q,HH,Hs) J ' dt J ds 

Jo(s, (jlj) U { j ] (Q£ - s, fx s , fij) S (Qt - £, fj, s ) ■ 

(24) 

In this scenario the jet and the soft function only receive 
contributions from the massless modes. The evolution 
factors employed in Eq. are just the ones from 
the massless theory. 



B. Scenario II: Q> M > QX> QX 2 

In this second scenario the mass modes enter SCET as 
long as their virtuality is above M. We therefore have 
to take them into account for the matching to the full 
theory at the hard scale /i# ~ Q and the RG-evolution 
above M. Thus the RG-evolution between /ih and /xm, 
which affects the current evolution only, accounts for the 
mass modes. When crossing the mass mode matching 
scale fiM ~ M we integrate out these massive modes. 
This leads to an additional current matching coefficient in 
comparison to the massless factorization formula. From 
there on everything proceeds as in the massless case. The 



factorization theorem thus takes the form 

-£ = Q 2 cr |C JJ (Q,M,/i H )| 2 £4 2) (Q,ti H ,ti M ) 

x \M H (Q,M,fi M )\ 2 U^ (Q,n M ,ns) Jdi 

J dsJ (s,fij)U^(Q£-s,fj,s,fij)So(QT-i,fj,s) . 

(25) 

In Eq. (|25|) the superscript (2) in the evolution kernel 
Utj (Q, hh, Hm) indicates that the evolution is now per- 
formed with the mass mode gauge bosons and the mass- 
less gluons. The hard matching coefficient acquires now 
a contribution arising from the non- vanishing SCET di- 
agrams involving virtual collinear and soft mass mode 
gauge bosons. The renormalized one-loop result of the 
corresponding sum of the virtual collinear and soft mass 
mode diagrams with on-shell external massless quarks 
reads 



5F eS (Q,M,n) = 




in agreement with Ref. [l4 |. The UV-divergences are 
mass-independent and agree with those from the mass- 
less diagrams. The calculation is carried out explicitly in 
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Sec. lVI Al and does not require any regularization in addi- 
tion to dimensional regularization. For the calculation of 
the mass mode collinear diagrams one has to account for 
soft mass mode bin subtractions. The result for the total 
renormalizcd Wilson coefficient at the scale fj, is given by 

C II (Q,M,») = C I (Q,M,»)-6F en (Q,M,(j,), (27) 

with C J (Q,M, n) being the hard current matching coef- 
ficient from scenario II in Eq. (f2~Tj). For M — > we now 

recover the correct massless limit, i.e. C 11 (Q, M, fi) 
2Co(Q,/i) with Co(<2>m) given in Eq. (0). The anoma- 
lous dimension of C II (Q, M, fi) is mass independent and 
just twice the one from the purely massless theory in Eq. 
©, so 

^ U< h {Q, »h,h) = 2 1H (Q, h)U ( h ] (Q, Mff j /*) • (28) 

The mass mode matching coefficient A4h(Q, M, Hm) 
is obtained by integrating out the massive modes at the 
scale ~ M. The result reads 

M H (Q, M, hm) = 1 + SF cS (Q, M, n M ) . (29) 

and is obtained by matching to the 0(a s ) full theory re- 
sult in Eq. (J5HJ for r < M 2 /Q 2 (M > QX), as explained 
in Sec. |Vl The matching correction in Eq. ((29)) agrees ex- 
actly with the sum of the collinear and soft mass mode 
corrections involved already for the hard matching con- 
dition of Eq. (|2"7| because the full theory results relevant 
for scenarios I and II are identical (being proportional 
to S(t)). This is related to the fact that the thresholds 
for mass mode production are located in scenarios III 
and IV so that the fixed-order full theory result (dis- 
cussed in detail in Sec. |Vj does not distinguish between 
scenarios I and II. In other words, the mass mode con- 
tributions SF e ff(Q,M, Hm), which are subtracted from 
the scenario I hard matching coefficient to obtain the 
infrared-safe C 11 (Q, M, fi) are added back in the mass 
mode matching coefficient when the mass modes are in- 
tegrated out. 

This is exactly the situation realized for hard coeffi- 
cients in the ACOT scheme P, 0- Since for [iu = P>H ~ 
Q, where the evolution factor U H (Q, /ih, Hh) is unity, 
the very same term is swapped between these two struc- 
tures, the continuity of the spectrum at the transition be- 
tween scenarios I and II is ensured up to higher order con- 
tributions which do not contain large logs and depend on 



the implementation of the fixed order terms. This shows 
that the transition between scenarios I and II has to be 
carried out for M ~ Q. Note that apart from the single 
cusp logarithm ln(— Q 2 /M 2 ), which is ensured to arise 
only as a single power to any order of perturbation the- 
ory, the mass mode matching correction 5F c ff(Q, M, /im) 
does not involve large logarithms for fj,M ~ M. So it has 
the correct properties of a mass mode matching condition 
being generated when the mass modes are integrated out 
for fi = Hm ~ M. On the other hand, the same terms 
generate the correct large infrared sensitive logarithms 
for fi = ~ Q to render the Wilson coefficient C 11 
infrared-finite in the limit M — > 0. 

Interestingly, the continuity of the transition from sce- 
nario I to scenario II is related to the fact that there are 
no higher order power corrections missed in the factoriza- 
tion formula of scenario II in the transition region M ~ 
Q. In this region the effective field theory contributions 
in the current matching coefficient C 11 (Q, M, fijj) and 
in the mass mode matching coefficient A4h {Q, M, /im) 
are proportional to SF e g(Q, M, hh) — 5F e g(Q, M, ^m) ~ 
{p-H - Hm)/Q ~ (Q - M)/Q < 1. This feature will re- 
main true at any order of perturbation theory in a s , but 
also concerning higher order M/Q power corrections in 
a subleading effective theory description for scenario II. 



C. Scenario III: Q > QX > M > QX 2 

The mass is between the jet and the ultrasoft scale. As 
indicated in Fig. [SI massive and massless collinear modes 
are part of the same collinear sector since the mass scale 
M is below their typical invariant mass QX. Thus mas- 
sive and massless collinear modes both contribute to the 
jet function. Between the jet scale /ij and the mass mode 
matching scale [Im the jet function evolves together with 
massive and massless collinear modes, and at \im the 
massive collinear modes are integrated out leading to an 
additional collinear matching coefficient which involves 
distributions. Concerning the current, its evolution and 
matching with respect to the mass mode matching scale 
Hm agrees exactly with the one of scenario II. Below the 
mass mode matching scale the current and jet evolution 
as well as the soft function agree with the ones obtained 
in the well known purely massless SCET setup. The form 
of the factorization theorem now reads 



da 
d7 



=Q 2 <j \C 11 (Q,M,hh)\~U ( h> (Q )\Mh(Q,M,» m )\ U£> (Q 

x JdlJdsJ ds'j ds" Jo +m (s,M lf ij)U { j 2) {s' - s,Hm,Hj)Mj(s" - s',M,^ M ) 
x U^\QI-s",iis^m)S {Qt-1^ s ) . 



2 rr(!) 



(30) 
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The matching coefficients C II (Q,M,fi) and 
■Mh(Q, M, hm) are the same as in scenario II, see 
Eqs. (f27| and ([29 ]) . 

The jet function Jo+ m (s, M, (i) contains additional 
contributions coming from the collinear mass modes (in- 
cluding soft mass mode bin subtractions) and has the 
form 



J 0+m (s,M, n) 



J (s,/i) + ^ irt (s,M,/i) 
+5J™\s,M). (31) 



The term <5 J™' (s , M , /i) contains only distributions and 
corresponds to virtual corrections and its renormalized 
expressions at one-loop order reads (s = s/ /i 2 ) 



-6 In 



a s C F 

47T 



9 - 2^ 



5(s) 



-41n J 



Sin 









s 



(32) 



The term 5J^ al (s, M) in Eq. (|33| contributes only when 
the jet invariant mass is above the threshold M and thus 
corresponds to real production of the massive gauge bo- 
son. At 0(a s ) it has the form 



sj;r\ s ,M) 



a s CF 
Air 



6(s - M 2 ) 



s \M 2 ) 



[2(M 2 - s)(3s + M 2 ) 



(33) 



Due to its physical character it is UV-finite and does 
not contain any explicit logarithmic /i-dependence. Fur- 
thermore, 5Jj° &1 (s,M) is zero at the threshold, so that 
no discontinuity arises due to real radiation. This fea- 
ture remains true also in the 0(a 2 ) massive quark cor- 
rection as the dispersion integration does not offset this 
property. For M — > the correct massless limit is 
reached by combining the real radiation pieces and the 
virtual contributions properly into distributions yielding 

Jo+m{s, M, fi) 2 Jo(s, /i). Details on the calculation 
of 5J™ Tt and 5J™ al are presented in Sec. IVIB1 Includ- 
ing the soft mass-mode bin subtraction the loop inte- 
grals are successfully regularized by dimensional regular- 
ization, and there is no need to introduce an additional 
regulator. The anomalous dimension coming from the 
virtual mass mode contribution 5 J?™ 1 is mass indepen- 
dent and exactly coincides with the one of the massless 
jet function in Eq. (fT4"|) . This yields 

d f 

^^C/j 2) (s,A*,Mo) = J ds'2~jj(s - s' , (i)Uy> (s' , (i, p, ) 

for the RG-equation of the jet function above the mass 
scale M, where the RHS is just twice the result obtained 
when the mass modes do not contribute, see Eq. ([3]). The 



(2) 

RG-equation for the current evolution factor U y H ' above 
M remains unchanged with respect to scenario II. 

The collinear mass mode matching coefficient 
M.j(s,M,hm) is obtained when the mass modes are 
integrated out from the jet function in its evolution 
down to the ultrasoft scale. At one-loop it reads 



Mj(s, M, pit) = 5(s) - W; irt (s, M, n M ) 



(35) 



and encodes just the virtual contributions contained in 
the mass mode jet function. Equation (|35[) is obtained by 
matching to the full theory result in the r region relevant 
for scenario III, M 2 /Q 2 < r < M/Q (QA > M > QX 2 ), 
as shown in Eq. (|57[) . In the collinear mass mode match- 
ing coefficient only terms corresponding to the distribu- 
tive virtual part of the mass mode jet function appear 
as required by consistency, since matching coefficients do 
not contain kinematic terms involving thresholds. For 
the thrust distribution the transition from the factor- 
ization theorem of scenario II to the one of scenario III 
should be made for values of thrust for which the jet scale 
is somewhat below the threshold of collinear massive real 
radiation, i.e. r < M 2 /Q 2 . 

Compared to the factorization theorem of scenario II 
we find that apart from the differences in the RG- 
evolution the distributive collinear mass mode contribu- 
tions in J™ t (s,M,y) are just swapped between the jet 
function and the collinear mass mode matching coeffi- 

(2) 

cient. Thus for jum = Hj where Uj(t — s" , hm, H>j) — 
5(t—s") the scenario III factorization theorem in Eq. ([50)1 
agrees with the one of scenario II in Eq. (|2"5|) (up to 
higher order terms which do not contain any large log- 
arithms and depend on the implementation of the fixed 
order terms) . This ensures the continuity of the transi- 
tion in the theoretical descriptions between scenarios II 
and III. We note that scenario III is based on the condi- 
tion M < Q\fr and no hierachy is assumed otherwise be- 
tween M and the jet scale. Thus the full M dependence is 
accounted for and scenario III has the same generic pre- 
cision for all allowed M values including the transition 
region to scenario II. 



D. 



Scenario IV Q > Q\ > Q\ 2 > M 



If the mass scale is below the ultrasoft scale the mas- 
sive modes are not integrated out at all because we have 
fj, ~ [is > M and the invariant mass of the ultrasoft fluc- 
tuations of order QX 2 exceeds the mass-mode scale M. In 
this scenario the typical invariant masses of the collinear 
and soft mass modes agree with the invariant masses of 
the collinear and ultrasoft massless degrees of freedom, 
respectively, and the mass scale M becomes merely a pa- 
rameter within the theory that is not relevant any more 
for the separation of modes. Therefore the collinear as 
well as the soft mass modes are contained in the mass- 
less collinear and ultrasoft sectors, respectively. In this 
scenario the massless and the massive degrees of free- 
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dom both contribute to all renormalization group evolu- 
tion factors as well as to the jet and the soft functions, 
and the factorization theorem does not contain any mass 
mode matching condition. The factorization theorem of 
scenario IV has the form 



Q 2 ao \C U (Q, M, hh)\ 2 U (2) (Q, hh, fi s 



da 
~dr~ 

x J d£ J dsJ 0+m (s,M,fj,j) 

X Uf\<Q£ - S, fls, fJtj) So+m {Qr - £, M, fig) ■ 

(36) 

Compared to the factorization theorems of scenarios II 
and III the hard current and jet mass mode matching 
coefficients have disappeared. Instead there is a new term 
in the thrust soft function due to the soft mass mode 
gauge bosons, 

+5S™\£,M). (37) 

In analogy to the mass mode contributions to the jet 
function the term SS™*^, M, /z) contains only distribu- 
tions and corresponds to virtual corrections. At one-loop 
order it reads (£ = £/fi) 




The term 8S!^(£, M) describes real radiation of the soft 
massive gauge bosons and has the form 



a s Cp 
4tt 



9(£- M) 



In 



M 2 



(39) 



It is UV-finite and does not contain any explicit in- 
dependence. We also find that 5S™ (£, M) is zero at 
threshold, which remains also true for the corrections in- 
volving secondary massive quark pair production. For 
M — » the correct massless limit is reached by combin- 
ing the real radiation pieces and the virtual contributions 

properly into distributions yielding So+m (A M, n) 
2So(£,fi). No mass singularities occur any more at this 
point. The anomalous dimension arising from the mas- 
sive contribution to the soft function is mass independent 
and coincides with the one of the massless contribution. 
This yields 

Hj-Ug\i,n,iMo) = J &l'2 ls {l-£',n)Tjf\£',n,n Q ) 

" (40) 
for the RG-evolution for the soft function above the mass 
scale M, where the RHS is just twice the result obtained 
when the mass modes do not contribute, see Eq. ([5]). 



However, since we have adopted the scale setting \x — fis 
here, no evolution of the soft function has to be accounted 
for. Notice that the evolution factors U H and U) in 



E. (|36[) remain unchanged with respect to scenarios II 
and III. 

The transition between scenarios III and IV is carried 
out for M ~ Qt. It is continuous in close analogy to 
the transitions between the scenarios II and III discussed 
above. However, the check is somewhat more involved 
here as it is based on the consistency relation 

2 Re [SF eS (Q, M, //)] 5(t) - Q 2 5J™\Q 2 t, M, 

-QSS^(Qr,M, t i) = (41) 

of the different virtual mass mode contributions. The 
relation can be verified by explicit analytic calculations, 
see the corresponding results given in Eqs. ([25)1 , ([521 and 
([55)) . On the other hand, we emphasize that the consis- 
tency relation ([41]) is part of an extended set of more 
general consistency relations that is associated to the 
freedom in fixing the renormalization scale \x to in prin- 
ciple any scale below Q. Relation (|4ip is also interesting 
by itself, since it expresses the fact that virtual mass 
mode contributions to the hard current, the jet function 
and the soft function are related to each other in a non- 
trivial way. We note that in practice the transition from 
scenario III to IV should be made for values of thrust for 
which the soft scale is somewhat below the threshold of 
soft massive real radiation, i.e. r < M/Q, We note that 
scenario IV is based on the condition M < Qt and no 
hierachy is assumed otherwise between M and the soft 
scale. Thus the full M dependence is accounted for and 
scenario IV has the same generic precision for all allowed 
M values including the transition region to scenario III. 

The way how the virtual mass mode contributions ap- 
pear and are interpreted depends on the choice of /i, 
as Eq. I|4ip enlightens. For example if we set /i = 
("bottom- up" approach, where the soft and the jet func- 
tions have to be evolved to hh), there are mass mode 
"matching" contributions for the soft and the jet func- 
tion, but none associated to the hard current contribu- 
tion. A more detailed discussion on the resulting consis- 
tency relations will be given in Ref. Q • 



V. FIXED ORDER FULL THEORY RESULT 

The mass mode matching coefficients given in the sce- 
narios outlined in Sec. IIVI are derived by matching the 
differential cross section in the different effective theory 
scenarios to the fixed order result ([56)) in the full theory. 
In this section we present the 0(a s ) full theory result 
by calculating the diagrams with massive gauge bosons 
shown in Fig. [3J We also discuss the various expansions 
relevant for the matching calculations in the different 
field theory scenarios. Since the full theory result has 
non-trivial threshold contributions related to the mass 
scale, these expansions deserve a separate discussion. 
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A. General Result 

The virtual contributions are known [23J and yield 
5F m (x) in Eq. (|2"2"j) using on-shell renormalization for the 
massless external quarks. Gauge invariance together with 
the on-shell normalization condition generate automati- 
cally SF m = at (p + p') 2 = [23]. What is left to be 
done, is the calculation of the diagram for real radiation 
of the massive gauge bosons and the integration over the 
three particle phase space. It is convenient to introduce 
the usual energy fraction variables, 



Xl 



2 Po 



X 2 



Q : 



£3 



2go 
Q ' 



M 
~Q 



where 



Xl 



x 2 +x 3 = 2 . 



(42) 



(43) 



Here p (p') is the quark (antiquark) momentum, M and 
q are the gauge boson mass and momentum, respectively, 
see Fig. [3] The double differential cross section in d = 4 
dimensions then reads 



1 dcr rcal a s C f 



1 



xl 



aodx\dx2 27r (1 — £i)(l — X2) 
2 (\ - 3x 2 + 2x^ 2 \ 4 

- X Z +I 4 + HH X 2 ) 

V l-Xi J 

The relation of the thrust variable r to £1,2,3 reads 

It -Pi I 

t = 1 — max 

t Q 



(44) 



= 1 — max ( xi,x 2 , \/ x 2 , — Ax 2 



(45) 



The resulting thrust distribution (including also the vir- 
tual contributions) displays structures to thresholds as- 
sociated to collincar and soft gauge boson radiation. It 
has the form 



1 da (x s Cf 
do dr 4-7T 



6F m (x)6(T) 



' mm I v \ ' max 



(a(t,x) + B(t,x) 
+ 9(t-t) 0(r max - r) C(r, x) \ , (46) 



where the minimal and maximal thrust values and the in- 
termediate threshold for the real radiation contributions 
are given by 



7m ax 

r 



1 

3 

x 



-l + 2y/l + ix 2 



(47) 
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FIG. 6: The phase space for the real radiation of a mas- 
sive gauge boson in terms of the variables x\ and X2 for 
x = M/Q = 0.2. The kinematic constraints yield the al- 
lowed area within the red, continuous lines. The blue dotted 
lines bound the areas with the same thrust relation to 11 and 
X2 and are given by x\ = xi, 2xi(2 — x-z) = (2 — X2) 2 — 4x 2 



and 22:2(2 



Xl) 



X2 — \ /f xT— 4a;' = 1 — r r 
Ei and E2 have minimal thrust configurations. At E3 the 
produced giuon is at res t with r = r. The triangle E1DE2 
corresponds to r = 1 — — Ax 1 and generates the function 
A(t, x). B(t,x) is obtained by increasing the areas of the 
segments E1E3D (E3E2D) to E1E4D (E5E2D) and integrat- 
ing with r = 1 — xi (t = 1 — X2)- Subtracting E5E4D with 
the corresponding opposite thrust prescription gives the soft 
function piece C(t, x). 



= (2 

4x 2 = 1 



x\Y — Ax' L . D (located at x\ = 
lx ) is the point of maximal thrust, 



The real corrections associated to the threshold at 
Tmin = M 2 / Q 2 correspond to collinear radiation, and 
those associated to the threshold at f = Mj Q correspond 
to soft radiation. As illustrated in Fig. |6l we have dis- 
tributed the contributions from the various phase space 
regions such that collinear and soft terms both extend to 
the endpoint at r max . This facilitates the identification 
of the singular terms necessary for the matching calcu- 
lations carried out in the SCET framework since in this 
way they are compatible with the contributions coming 
from the jet and the soft functions. Note that in Eq. (|46"]) 
we have factored out an overall loop factor and rede- 
fined a s CF/2Tr5F m (x) = Re[5F m (x)] with SF m (x) given 
in Eq. The functions A(t,x), B(t,x) and C(t,x) 

are obtained by integrating da rcal /dxi dx 2 in the phase- 
space regions as described in the caption of Fig. [5] The 
analytic results read 
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A(t, x) 



B(t,x) 



C(t,x) 



>-; TKl ~ r) [1 - w 2 2r(l + 2w) + 5r 2 ] 
tw{w — T) J 

+ — ^ \9 - 12w + Uw 2 - Aw 3 + w 4 + 4t(3 - 2w + w 2 ) - 2r 2 (l - w) 2 - At 3 + t 4 1 In 
2u> 2 L J 

(1 + w - 3r)(l + r - w) 



(48) 



5t 3 (w — t) 



1 



+ _L[ (3 + ^-2t 2 (1 



[tu(l - u) 2 ) - t(1 + 14w - bw 2 + Aw 3 ) + r 2 (7 - Aw){2 - w) + r 3 (3 + Aw) - At 4 
At{w — t) 



r 4 lln 



w 2 — (1 — 



— [1 - w 2 - 2r + 5r 2 l [l - w 2 - 2r(7 + 2w 2 ) - llr 2 + 4r 3 l 
8r J J 



l r 

27 



(3 + tI; 2 ) -2r 2 (l + ^ 2 ) + r 



In 



4r 2 



w 2 — (— 1 + t) 2 
I 



(49) 



(50) 



with 



4a; 2 . 



(51) 



The result for fixed-order thrust distribution in arbitrary 
units is displayed for x = 0.2 in Fig. [7] As a cross-check, 
from Eq. PS)) we obtain the correct massless expression 
for x -> 0, namely [M [H| 



1 da real m^o a 8 C f 



(TO 



dr An 
~6{t) lnr" 











r 



1 - r 



[6 + 3r - 9r 2 



- - 6 + 6t ) ln(l - 2r) 



(2-4r) lnr 
(52) 



valid for < r < 1/3. To obtain this result involving the 
plus-distributions it is required to exactly account for the 
analytic behaviour of the collinear and soft thresholds 
when taking the massless limit. 




0.2 0.3 0.4 0.5 



FIG. 7: The three contributions A(t,x), B(t,x) and C(t,x) 
to the 0(a 3 ) full theory thrust distribution for the radiation 
of a massive gauge boson for x — M/Q = 0.2. The sum 
of all terms, given by the thick green line, is positive. The 
vertical green line at r = indicates the virtual corrections 
proportional to <5(t). 



B. Expansions 

The massless fixed order thrust distribution shown in 
Eq. (|52j) only depends on r, and the expansion param- 
eter relevant for approaching the dijet limit is simply 
A ~ \fr ^ 1. There is only one single threshold located 
at t = 0. Thus in the massless limit we can identify the 
singular pieces that can be used for the matching calcula- 
tions within SCET by performing an expansion in small 
r, and the collinear and soft contributions are not sepa- 
rated. In the presence of the massive gauge bosons the 
collinear and soft thresholds are separated due to their 
different mass-dependence: collinear radiation starts at 
Tmin = M 2 /Q 2 and soft radiation starts at f = M/Q. 
Thus the identification of the singular terms - which shall 
include the full mass-dependence - needs to account for 
the location of the thresholds. In this subsection we dis- 
cuss the required expansions for the full theory results 



and show a simple example how it is applied in a match- 
ing calculation. 

The S(t) term is left unexpanded since we aim at de- 
termining the full mass-dependence of the singular con- 
tributions. The term proportional to 9{t — x 2 ) contains 
the collinear contributions. Here the jet invariant mass 
Qyfr and the mass M have to be considered at the same 
footing. Since the mass threshold is treated exactly, we 



expand for r 



<C 1. Here only the function B con- 



tains singular terms as it originates from the collinear 
regions in the X\-X2 phase space where ii ~ 1 or i 2 ~ 1: 



B{t,x)^* 



(x 2 -T)(3r + x 2 )+4r 2 ln(^) 
+ O(t ,x°). 



(53) 



The function B and all terms in the expansion vanish 
at t = x 2 giving a continuous turn-on of the singular 
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C(t _1 ) terms (with r ~ x 2 ). On the other hand, the 
function A{t,x) does not contain singular C(r _1 ) terms 
as it stems from configurations which are neither collinear 
nor soft. Interestingly the function A is non- vanishing at 
and responsible for the step visible in Fig. [7] at the 
c 2 <C 1 the expansion for the function 



T = X 
T = 

A reads 



A(t,x) t - 



rln r 



+ 0(t,x 2 ) 



(54) 



showing that it contributes at O(t ~ 1). Although the 
contribution from function A exceeds the singular contri- 
bution from B numerically, it only contains non-singular 
terms which are not described in the factorization the- 
orems. It remains to discuss the term proportional to 
6{t — x), which contains the soft contributions. Here the 
soft scale Qt and the mass M have to be considered at 
the same footing. Since the soft mass threshold is treated 
exactly we expand in r ~ x <C 1 • For the function C this 
yields singular C(r _1 ) terms which read 

C(t,x) ^ -— In(-) +O{t ,x°). (55) 

T \X / 

Similar to B, the singular 0(1 /r) terms turn on contin- 
uously at r = x. 

Assembling all singular pieces to be considered for the 
matching calculations we obtain 



1 da 



full th. 



Co dr 
a s C F 

47T 



8{T)+2Re[SF m {M/Q)} S(r) 



FO 



9(Q 2 t - M 2 



2(M — Q t)(3Q t + M 



f2\ 



— In. 

r V M 2 



r -M)-ln( 

t V M 2 



(56) 



The terms corresponding to the collinear and soft thresh- 
olds yield functions which match exactly to the real ra- 
diation effective theory contributions from collinear (see 
<5J™ al in Eq. ([33])) and soft (see 5S*™ 1 in Eq. figft) mass 
modes, respectively. We note that from Eq. (|56|) the sin- 
gular terms in the massless limit (first three terms on the 
RHS of Eq. (1521) ) can be recovered, requiring the proper 
combination of the terms into plus-distributions and 5- 
functions in r. 

An interesting issue we would like to mention is that 
the parametric precision of the expansions leading to the 
singular contributions in Eq. (j5"u| is not uniform and de- 
pends on t. Above the collinear threshold (t > T m - U1 ) 
the expansion is valid up to higher order power correc- 
tions of order M 2 /Q 2 . On the other hand, above the 
soft threshold (r > f ) the expansion is valid up to higher 
order power corrections of order M/Q. This is an in- 
trinsic property and also inherited to the factorization 
theorems we derived in the various scenarios. So, from a 
strict power counting point of view, the soft sector would 
need to be treated with the first order power corrections 
included to reach the same parametric precision as the 



collinear sector. We stress, however, that this does not 
affect in any way the consistency of treating only the sin- 
gular collinear and soft mass mode contribution in the 
effective theory description. In practice this might make 
the treatment of subleading power corrections more im- 
portant for the soft mass modes than for the collinear 
ones. 

As an example for the matching procedure from the 
underlying theory to SCET we now derive the one-loop 
collinear mass mode matching contribution JV^p in sce- 
nario III where M 2 /Q 2 < t < M/Q. Expanded in fixed 
order (i.e. setting /i = /in = U>J = = Ms); the mass 
mode contributions to the factorization theorem for sce- 
nario III in Eq. (|30|) read 



1 da SCET 



CO 



dr 



in 



= 5{t) +2 Re [SF m (M/Q)]S(T) 



+ Q A 6 J™ 1 (Q 2 t, M, ii) + M { P (Q z t, M, fi) 



V- M 2 )Q 2 5J™ a \Q 2 T, M) 



To determine A4 y we have to compare this result to the 
expanded full theory result of Eq. (|56| for M 2 /Q 2 < r < 
M/Q. This yields the condition 

M { }\Q 2 r, M, fi) = -SJ^(Q 2 t, M, fi) . (58) 

canceling exactly the distributive pieces of SJ m as de- 
scribed in Sec. IIV CI All other mass mode matching co- 
efficients are derived in a completely analogous way. 



(57) 



VI. CALCULATIONS WITH MASS MODES 

In this section we describe the effective theory calcula- 
tions for the mass mode contributions to the infrared-safe 
hard Wilson coefficient and the jet and the soft functions 
entering the factorization theorems in Sec. IIVI For the 
regularization of all Feynman diagrams we use dimen- 
sional regularization. The integration measure in light- 
cone coordinates reads 



d d k 
(2^p 



-d/2 



2 2vr 2vr 



\ 2 i 



Ak i k 



d-3 



(59) 



with k± = \kj_\ for integrals not depending on the angles 
between the transverse momenta. An important techni- 
cal point is that we encounter divergences in single di- 
agrams which are not regularized by dimensional regu- 
larization. In our calculations we do not employ an ad- 
ditional regulator (see e.g. Refs. because these 
divergences turn out to cancel among the proper set of 
contributions arising from diagrams together with the re- 
spective soft mass mode bin subtractions (required to 
avoid double counting). We stress again that the collinear 
and the soft mass modes have the same typical invari- 
ant mass, while for massless modes the invariant masses 
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FIG. 8: Non- vanishing EFT diagrams for the computation of 
the hard matching coefficient, soft mass mode bin subtrac- 
tions are implied for the collinear diagrams. 



for collinear modes exceed the ones for ultrasoft modes. 
In this respect the soft mass mode bin subtraction is 
somewhat different from the zero-bin subtraction [26l - [28| 
for massless collinear diagrams, which can be related to 
vanishing large collinear label momenta. We note that 
we have checked explicitly that the contributions com- 
ing from longitudinal polarizations of the massive gauge 
bosons vanish due to gauge invariance. To achieve this for 
the calculations involving collinear mass modes the inclu- 
sion of soft mass mode bin subtractions turned out to be 
crucial. So the soft mass mode bin subtractions are essen- 
tial to isolate gauge-invariant structures and to achieve 
factorization, as has also been stressed in Ref. [3| ■ In the 
following we present all calculations in Feynman gauge. 



A. Vertex Corrections and Hard Current Matching 

In this section we discuss the calculation of the 
collinear and soft mass mode contributions to the effec- 
tive theory vertex corrections feeding into the matching 
for the hard current Wilson coefficient at the matching 
scale fin ~ Q shown in Eq. (j26[). The result is valid 
for all scenarios where the mass scale M is below the 
hard scale Q (see Eqs. ([21]), ([3D]) and (|55]l). This can 
be understood from the fact that the way how collinear 
and soft mass modes are distributed with respect to the 
massless modes below the hard scale Q does not affect 
the hard contribution itself. 3 We note that the same re- 
sult has been obtained before in Refs. [l3| and [l4| using 
an additional analytic regulator and the A-regulator, re- 
spectively Here we show the calculation once more using 
only dimensional regularization and the fact that all sin- 
gularities cancel once soft mass mode bin subtractions are 
included. We also clarify the role of the soft mass mode 



3 We carry out the calculations for external massless collinear 
quarks which are on-shell. This leads to exactly the same calcu- 
lations for all scenarios with M < Q. 



bin subtractions for the collinear wave- function renormal- 
ization factors. The relevant collinear and soft diagrams 
with massive gauge bosons to be calculated are shown 
in Fig. [5J where p (p') denotes the momentum of the 
collinear massless external quark (antiquark). 

Applying the SCET Feynman rules gives for the soft 
mass mode diagram 



V s =2ig 2 jl 2e C F 
P 



d d k 



P 



(2w) d [(p - fc) 2 + ie] 
■'+ -J- fc+ 1 
" [(p> + fc) 2 + ie] k 2 - M 2 + ts ^ ' (60) 

Using p — (0, Q, 0) and p' = (Q, 0, 0) and the soft scaling 
k^ ~ Q(Am)AmjAm) for the soft momentum we obtain 
for Am <. 1 



V s = 2ig 2 jl 2£ C F 



d d k 



1 



1 



(2it) d [-k+ + ie] [k~ + ie] 



[k 2 — M 2 + 



ie 



(61) 



The expression in Eq. (|fJT|) can as well be obtained di- 
rectly from the SCET Feynman rules upon field redefini- 
tion such that the soft mass mode gauge bosons couple to 
the collinear massless quarks through mass mode Wilson 
lines 0, [HI ■ The n-collinear diagram yields 



V n = 2ig 2 C F fi 2e 



d d k 



{2ir) d [{p - fc) 2 + ie] [fc- + ie] 

1 



[fc 2 - M 2 + ie] 



and using the form of the external momenta (and likewise 
the collinear scaling ~ Q(X 2 4 , 1, Xm)) we obtain 



V n - 2ig 2 C F fi 2e 



d d k 



Q 



(2ir) d [k+{k- -Q)-k\ +ie] 



[fc- +ie][k 2 - M 2 + ie 

2ig 2 C F p, 2<i I n £ n ,pF(;n,p' ■ 



(63) 



The n-collinear diagram V n is obtained in an analogous 
manner and involves the function I n ■ It is obtained from 
the n-collinear diagram by swapping the plus- with the 
minus-components as well as p with p' . Expanding the 
?i-collinear diagram of Eq. (1521) in the soft mass mode 
regime fc ~ Q{\m, Am, Am) we obtain the soft mass mode 
bin subtraction contribution, 



V, 



n.OM 



2ig z C F fi 



,~,2t 



d d k 



Q 



(2ir) d [-Qk+ + is] 



ie\ 



1 



[fc 2 - M 2 + ie] 

2ig 2 C F fl 2e I n ,0M £n,pFt;n,p 



Cn.p^Cn,p' 



(64) 
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with the analogous expression for the n-collinear sector, 
which we call I n ,oM- Comparing with Eq. (|0 1 [) we find 
the known relation I n flM — In.oM — h and the sum of 
all vertex diagrams including the soft mass mode sub- 
tractions gives |14| 

SF eS ~ (In-In,QM) + (In-Ifl,QM)+Is = J n +Ift-J a (65) 

as for the massless case with 0-bin subtractions [13] ■ 

We compute the integrals /„ , I ft and I s by first carrying 
out the &; + -integration with the method of residues and 
then the ^-integration. So the integrals /„ and J fi are 
not treated symmetrically and lead to different analytic 
expressions for the remaining k~ integral. In the fc~ in- 
tegration all singularities cancel properly in the sum of 
all terms and the result is finite in dimensional regular- 
ization. We start from the integral /„ given in Eq. (1551) . 
Carrying out the k + - and fcj_-integration as described 
above we obtain 



(66) 

with z = k /Q. The integral in Eq. (|66[) is divergent 
showing that dimensional regularization fails. For the /„ 
integral it is convenient to split it up as 



In .=/« - fr 2e B (o,M 2 ,o), 



(67) 



where Bo(p 2 ,mf,m|) denotes the usual two-point func- 
tion. After performing the k + - and fc^-integrations we 
arrive at (x 2 = M 2 /(Q 2 + iO)) 



l-oo 1 

< / dz — 
/o 



d/2-2 



(68) 



The result differs from Eq. (|66)) since the integrations for 
I n and I n were performed in an asymmetric way. Finally 
we compute I s given in Eq. (I5T|) in the same way and 
obtain 



(47r)d/2 n 2 "2)( M2 ) d/2 2 



dz 



(69) 



which is again not regularized by dimensional regulariza- 
tion. 



Upon summing all contributions according to Eq. (|6"5)) 
the singularities cancel exactly. For this purpose we split 
I s given in Eq. (|69p into an integration from to 1 and 
from 1 to oo. The first integral is then recombined with 

For the sum of the divergent 



integrals we then obtain 
- 1 (1 - z) 11 ' 2 - 1 - 1 



z 

l-co X 2 + Z 
'1 

H d 



dz + 



d/2-2 



1 1 



d/2-2 



-dz 



(x 2 + z)z 

ri'fWi-fl-W-." 



(70) 



where H a denotes the Harmonic number function. 

Finally, we compute the wave function renormaliza- 
tion diagrams Wf >a and Wf t b for the external massless 
collinear quarks. The sum of both self-energy diagrams 
for arbitrary collinear external momentum p^ can be 
readily combined and reads 

W f , a +W f , b = -g 2 C F (d~2)*h 2e ^ f ^ 



2 ^ Q J (2tt 



P 2 + Qk^ 



[p 2 (l + k-/Q) + Qk+ + k 2 + ie][k 2 - M 2 + ie] ' 

(71) 

which is known to reproduce full theory self energy graph. 
The wave-function renormalization contribution can be 
identified in the limit p 2 —> (and e — > 0) giving 

.ff, p 2 



Wf, a + W f , b 



i- — W f 

2p- 1 



where 



W f 



■in 



(72) 



(73) 



The soft mass mode bin subtraction for the wave func- 
tion contribution turns out to be power suppressed by 
M/Q. To demonstrate this let us first consider the 
power counting of the collinear mass mode self-energy 
diagrams given in Eq. ()T1 [) : the collinear loop momenta 
scale like k^ = (k^,k + ,k±_) - (Q,M 2 /Q,M) and the 
external massless collinear modes interacting with the 
collinear mass modes thus obey the same scaling giving 
p 2 - (p + k) 2 - M 2 . This yields Wf, a + W f , b - M 2 /Q 
for the result in Eq. (JTTJ) and gives the same count- 
ing for the wave- function contribution given in Eq. (|72|). 
p 2 /Qx W f - M 2 /Q. This demonstrates that Wf ~ 0(1) 
as we have obtained also in the explicit result shown in 
Eq. (l73|) . For the mass mode bin subtraction to Eq. ([7Tj) 
we have to apply the counting fcg ~ (M, M, M) for the 
loop momentum and this increases the typical invariant 
mass of the massless collinear quarks that interact with 
the soft mass mode gauge bosons top 2 ~ (p+k s ) 2 ~ QM . 
The mass mode bin integral that emerges therefore has 
the form of a simple tadpole graph, 

{Wf, a + W /)6 ) m = -9 2 C F (d - 2)*f M 2e ^ 



/„, the second with /, 



d d k 



1 



(2ir) d [k 2 - M 2 + it 



(74) 



17 



Ja 



® ® 




® 



FIG. 9: Non- vanishing EFT diagrams for the computation of 
the jet function. The required soft mass mode bin subtrac- 
tions are implicit. Concerning J a also the right-symmetric 
diagram has to be taken into account. 



with the counting (W>, + W ftb ) M ~ M 2 /Q. The 
same counting also applies to the resulting wave-function 
renormalization contribution, i.e. p 2 /Q X (Wf)oM ~ 
M 2 /Q, and because p 2 ~ QM we find that (W/)om ~ 
M/Q. This shows that the mass mode bin subtraction 
for the collinear wave function renormalization due to the 
mass modes belongs to a subleading treatment beyond 
the scope of the treatment discussed here. 4 Combining 
Eqs. ([q7)1 . ([70)1 and ([75)1 and expanding in e we arrive at 
the final result 



SF^ e (Q,m,fx) = 
=2ig 2 C F M 2£ 

47T 



2 In 



-0 



2 



L 



'M 2 



S (0,M 2 ,0) 
M 2 



In 



In 



/' 

57T 2 

6 



(75) 



In Eq. (175)1 the UV-divergences do not depend on the 
mass and agree exactly with those of the sum of the cor- 
responding collinear and soft effective theory diagrams 
for massless gluons, see Eq. ©. After renormalization 
we thus obtain the result for 5F c g given in Eq. (|26|) . 



The mass mode contributions to the jet function arise 
in the situation A > Am, where the collinear mass 
modes have the typical invariant mass QX of the mass- 
less collinear modes. So we can use the scaling k^ ~ 
Q(X 2 , 1, A) for the n-collinear mass mode gauge bosons. 5 
The jet momentum and its invariant mass are denoted 
with p^ and s = p 2 = p~p + = Qp + , the prescription 
s — > s + iO is understood. We keep p + explicitely non- 
zero in the subsequent computations to account for the 
non-vanishing off-shellness s of the jet sector and use a 
frame with p± = 0. 

Taking into account the traces and the prefactor of 
the jet function matrix element in Eq. we obtain for 
diagram J a 



2g 2 C F _ 2e f d d k p- -k 

J a — M 



(27r) d [(p - k) 2 + ie] 
1 



[k~ + ie][k 2 - M 2 + is] 



(76) 



For our choice of the external momenta this yields 



J a = 



2g 2 C F _ 2( 



TVS 



d d k 



Q-k- 



(2ir) d [s(l - k-/Q) - Qk+ 
1 



■ ie\ 



[k~ + ie][k 2 



M 2 



■ ie\ 



2g 2 C F _ 



(77) 



The corresponding soft mass mode bin subtraction with 
the scaling k 11 ~ Q{Xm,Xm, Am) for scenario III (or k 11 ~ 
Q(X 2 , A 2 , A 2 ) for scenario IV) gives 



B. Jet function 

In this section we calculate the 0(a s ) mass mode gauge 
boson contribution to the n-collinear massless quark jet 
function J m . n (s, M, fi) and relate it to the thrust jet 
function of Eq. (1311) appearing in the factorization the- 
orems (|3U)1 and (|3"6")1 for scenarios III and IV where 
A > Am- The diagrams contributing to the n-collinear 
jet function are listed in Fig. |H1 Similarly to the hard 
function calculation, the singularities that arise in the 
collinear mass mode contributions and which are not han- 
dled by dimensional regularization are canceled by the 
soft mass mode bin subtractions. 



2g 2 C F „ 2t f d d k Q 

Q '° M its ^ J {2Tr) d [s-Qk+ +ie}[k~ +ie] 

1 

X [k 2 - M 2 + ie] 

= H Ja,0M- U°) 

ITS 

Note that for scenario III the invariant mass term s is 
included as an infrared scale to account for the non-zero 
jet-invariant mass. In close analogy to the wave func- 
tion contributions to the vertex corrections discussed in 
Sec. I VI Al the sum of the self-energy diagrams J& and J c 



4 Since M/Q-suppressed terms do not exist in the non-singular 
collinear terms that can be obtained in the full theory calcula- 
tion, the subleading effective field theory treatment contains a 
mechanism that makes these contributions vanish identically. 



J In scenario III the soft mass modes have typical momenta k ~ 
Q(Ajvfi AM: Am)- In scenario IV the soft mass modes have typical 
momenta k ~ Q(A 2 , A 2 , A 2 ). 
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reads 



The sum of all diagrams in the limit e — > yields 



Jh + J c — 



g 2 C F {d-2)„ 2t 
"A* 



us 



d d k 



1 



2n) d [fc 2 - M 2 + ie 
1 + Qk+/s 



[s(l + k-/Q) + Qk+ + k 2 +is] 



g 2 C F (d~2)„ 



TVS 



(h + h) ■ 



(79) 



As we have discussed for the vertex corrections in 
Sec. IVI A[ the soft mass mode bin subtraction terms to 
the self-energy diagrams belong to a subleading effective 
theory treatment and are thus not considered. 

To compute the integrals, we apply the same technique 
as in Sec. IVI Al So we first carry out the k + and k± in- 
tegrations and combine the contributions of all diagrams 
prior to the k~ integration. For I a we then arrive at 



— (4^ r l 2 -2, )(M2)d/2 1 



■ l" 1 — (1 — z) d / 2 ~ l ( I 
/o z 



~M 2 - Z ' 



(80) 



with z = k~/Q. There is an unregularized divergence 
for z — > 0, which will be cured by the corresponding soft 
mass mode bin subtraction, 



W=(j^r(2-|)(MV/" 

f°° dz ( S \d/2-2 
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(84) 



We can take the absorptive part with the help of the re- 
lations in the appendix of Ref. 6], which leads to the 
unrenormalized mass mode contributions to the jet func- 
tion 



8(a) 
4 In 



4tt 
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(85) 



We see that the UV-divergences are mass-independent 
and agree with those from the purely massless jet func- 
tion. Multiplying the result by a factor of 2 to account 
for the combination of the two hemisphere jet function 
into the thrust jet function and using the jet function 
renormalization counterterm given in Eq. (1131) we get 
SJ m (s,M,n) in Eq. 



Subtracting Eq. ([51) from Eq. flU]) yields 



S ^/2-2 



OO 7 

*,1 



S \d/2-2 



(82) 



which is finite in dimensional regularization. 

Next, we proceed to the calculation of the self energy 
integrals If, + I c . Defining z = —k + /p + we arrive at the 
expression 



x I dz 
lo 



('-^'('-ai') 



o sd/2-2 



C. Soft function 

In this section we calculate the 0(a s ) contributions of 
the soft mass mode gauge bosons to the soft function for 
scenario IV, where A 2 > Am and the mass mode scale 
is below the ultrasoft scale. The non-vanishing diagrams 
are displayed in Fig. [101 Diagram S a (Sb) corresponds 
to virtual (real) corrections. Their contributions to the 
thrust soft function read 



8{i) 



, (86) 



(fc+ + ie)(k- - ie)(k 2 - M 2 + ie) 
S b =4ng 2 C f ^J ^ 6(fc+ + k-)6(k 2 - M 2 ) 

Q(k~ - k+)5{£ - k+) + 6(fc+ - k-)S{£ - fc-) 



(fc+ + ie)(k~ - ie) 



(87) 
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FIG. 10: Non-vanishing Feynman diagrams for the computa- 
tion of the one-loop mass mode contributions to the thrust 
soft function. Corresponding symmetric configurations are 
implied. 



Although the sum of all contributions is finite in dimen- 
sional regularization, there are again individual singular 
regions with respect to large (or small) light cone mo- 
menta k + and k~ , or in other words with respect to ra- 
pidity. These singularities, however, cancel themselves 
and no manipulations or subtractions of any kind are re- 
quired. To perform these integrals it is convenient (but 
not mandatory) to introduce an additional IR regula- 
tor besides dimensional regularization. We choose a- 
regulators motivated by Ref. [2j|, which are additional 
analytic regularizations for the light cone k compo- 
nents, 



dk- 



dfc H 



dk~ 



\k- 



ffl 



In this context the scales vi% are simply auxiliary scales 
to maintain the dimension of the integral. We stress that 
these scales do not carry any physical meaning and are 
not related to any kind of renormalization as they drop 
out entirely from the final result for the soft function. 
Implementing these replacements, S a and Sb are indi- 
vidually well defined. The virtual diagrams S a eventu- 
ally yield scaleless integrals and give zero, similar to the 
case of a massless gluon exchange (see Ref. |6[). Using 
v\ = V2 = fx for convenience the real radiation diagram 
S d gives 



n ,,2e+ai+a2 p 7Ee 
Sh = Q l° f ^ — % M -2 ai -2e 9 ^ e -l+ ai - a2 



2tt r(l-e) 
T(l - £ )r(ai +e) 

r(i + «i) 

+ (oti <-> a 2 ) , 



+ B M 2 /e 2 (ai + e, 1 - e) 



(89) 



where B z (a, b) denotes the incomplete Beta function. Af- 
ter expanding in a% and then a 2 , taking the limit e — > 



and adding the symmetric diagram one obtains (i 
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for the thrust soft function. We see that the UV- 
divergences are mass- independent and agree exactly with 
the well known divergences of the massless gluon contri- 
butions to the soft function. Subtracting the divergences 
with the counterterm of Eq. (JTSJ) leads to the renor- 
malized expression of the thrust soft function given in 
Eq. (021). 



VII. SECONDARY MASSIVE FERMIONS 

As outlined in the introduction our main motivation 
for studying the effective field theory setup involving the 
mass mode gauge bosons is the treatment of massive 
secondary quark- antiquark pairs arising from the 0(a 2 s ) 
gluon splitting diagrams of Fig. Q] For cases where the in- 
variant mass of the quark pair enters the observable the 
corresponding effects can be calculated from the mass 
mode gauge boson results with the help of dispersion in- 
tegrations. In this section we outline the method. In 
our presentation we assume the existence of rif massless 
quark flavors and one quark species with mass m. The 
presentation of calculational details for thrust and other 
quantities is beyond the scope of this paper and will be 
given in subsequent publications. 

The vacuum polarization Tl(q 2 ) describing a massive 
fermion bubble can be expressed through an integral over 
its absorptive part. The unsubtracted (unrenormalized) 
version reads 



IT(<Z 2 



dM2 T72 2— 

4m 2 M 2 - q l - it 



Im [n(M 2 )] , (91) 



while the subtracted (on-shell and finite) version has the 
form 



n°V) = u( q 2 ) - n(o) 
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The absorptive part in d dimensions reads 
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The correction of the massive quark-antiquark loop to 
the gluon propagator, as illustrated in Fig. [3J can be 
expressed in terms of dispersion relations. Denoting the 
external momentum by the relations read 



KM 2 



(-i) 2 g, p W°(g 2 )g a „ 
(q 2 + ie) 2 

1 r dM 2 



Im [II(M 2 )] 



M 2 



(94) 



for the unsubtracted and unrenormalized version and 

n^ os (g 2 ) 
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n Jim 2 



M 2 



q 2 - M 2 + ie 



(95) 



for the subtracted version involving the on-shell renor- 
malized vacuum polarization. 

The relations in Eqs. and (pj)) involve massive 
gluon propagators 6 and provide the formal connection 
to the mass mode gauge boson considerations discussed 
in this work. To determine the corrections arising from 
the secondary massive quark-antiquark pair one can first 
calculate the corresponding 0(a s ) diagrams with mas- 
sive "gluon" propagators and then convolute the result 
according to Eqs. (|94[) and (|95|). We note that in general 
the final convolution needs to be carried out involving en- 
tirely d dimensional expressions for the absorptive part 
of the vacuum polarization as well as for the 0(a s ) di- 
agrams with the massive "gluons" . The various forms 
of the factorization theorems discussed in Sec. |TV] con- 
cerning the arrangement of matrix elements, matching 
conditions and renormalization group factors remain. 

Massive secondary quarks contribute to the renormal- 
ization evolution only in those terms of the factorization 
theorem dealing with scales above the matching scale 
Mm ~ T7i. The dispersion integration method can deal 
with this issue as well. Indeed, the way secondary quarks 
contribute to the renormalization group evolution can be 
related to the use of the subtracted versus the unsub- 
tracted dispersion relation: when the massive secondary 
quark pair represents an active flavor contributing to the 
renormalization group evolution together with the mass- 
less quarks and gluons, it is convenient to employ the 
unsubtracted dispersion relation of Eq. (|94)l . This is 
the case for the calculations of the massive secondary 
quark pair effects to the hard Wilson coefficient in sce- 
narios II, III and IV, to the jet function in scenarios III 
and IV and to the soft function in scenario IV. Here 



the massive quarks contribute as dynamical degrees of 
freedom and enter in the renormalization group evolu- 
tion of the respective components in the factorization 
theorem as well as of the strong coupling. In particu- 
lar we find that the associated UV- divergences are mass- 
independent and agree with the ones obtained for mass- 
less secondary quarks (as we also found in an analogous 
way for the one-loop corrections from the mass mode 
gauge bosons). The use of the unsubtracted dispersion 
relation also implements mass-dependent logarithms in 
the UV-finite terms that are essential to reach the correct 
massless limit for m — > 0. Thus in these cases the out- 
come concerning renormalization group evolution for one 
massive secondary quark-antiquark pair is simply that 
the number of active flavors is nt + 1, where nf is the 
number of massless quark flavors. In the factorization 
theorems discussed in Sec. lIVI this affects all renormaliza- 
tion group evolution factors which included massless as 
well as massive gauge boson contributions (indicated by 
the superscript (2)). All other renormalization group fac- 
tors (indicated by the superscript (1)) evolve with n/ ac- 
tive flavors, since they describe evolution for scales where 
the massive quarks have been integrated out. 

When the massive secondary quarks are integrated out, 
one has to employ the subtracted dispersion relation of 
Eq. (|95[) with the vacuum polarization being in the on- 
shell scheme, thus giving the correct decoupling limit. In 
this situation the massive quark flavor is not dynami- 
cal and does not contribute to the renormalization group 
evolution. This is for example the case for the calcula- 
tions of the massive secondary quark pair effects to the 
hard Wilson coefficient in scenario I. The difference be- 
tween using the subtracted and the unsubtracted disper- 
sion relation, which is related to the vacuum polarization 
function at zero momentum transfer, 11(0), correctly im- 
plements the matrix element corrections related to the 
matching relation that connect the strong coupling a s 
schemes with n/ and nt + 1 dynamical flavors, i.e. the 
well-known decoupling relation. We note that in the per- 
turbative expansion for the mass mode matching condi- 
tions for the current (Mh in the factorization theorems 
of scenarios II and III) and the jet function [M. j in the 
factorization theorem of scenario III) one can either use 
the unsubtracted or the subtracted versions of the dis- 
persion relations depending on whether one adopts the 
rif + 1 or the n f flavor running for the strong coupling in 
the perturbative series for the matching functions. The 
previously mentioned corrections related to the difference 
between using the subtracted and the unsubtracted dis- 
persion relation ensure that the prediction is independent 
of the scheme choice for the strong coupling. 

We conclude with a brief discussion of the numerical 
size of secondary bottom quark mass effects in the singu- 
lar partonic thrust distribution at N 3 LL order 7 where all 



The longitudinal contributions differ, but do not contribute due 
to gauge invariance. 



7 We use the canonical SCET counting, where N 3 LL order refers 
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FIG. 11: (a) The partonic N 3 LL order thrust distribution at 
Q = 14 GeV with secondary bottom quark mass effects for 
mj = 4.2 GeV (blue solid line). The curve describes primary 
production of nj = 4 massless quarks and secondary pro- 
duction of the massless quarks and the bottom quark. The 
red dashed line shows the result for vanishing bottom quark 
mass, (b) Relative size of secondary bottom quark mass cor- 
rections. The black dotted lines indicate the boundary of 
theory description in scenarios III and IV. 
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FIG. 12: (a) The partonic N 3 LL order thrust distribution 
at Q = 500 GeV with secondary top quark mass effects for 
m t — 175 GeV (blue solid line). The curve describes primary 
production of n/ =5 massless quarks and secondary produc- 
tion of the massless quarks and the top quark. The red dashed 
line shows the result for vanishing top quark mass, (b) Rel- 
ative size of secondary bottom quark mass corrections. The 
black dotted lines indicate the boundary of theory description 
in scenarios II, III and IV. 



two-loop mass corrections were calculated from the mas- 
sive gluon results via the dispersion method described 
above. 8 In Fig. [TTk the blue solid line shows the re- 
sult for the production of n/ = 4 massless primary 
quarks at Q — 14 GeV with secondary production of the 
rif = 4 massless quarks and the bottom quark with mass 
nib = 4.2 GeV. The red dashed line shows the result when 
the bottom mass is zero. We use profile functions similar 
to the ones used in Ref. [25| to describe the thrust depen- 
dence of the jet and soft scale. The massive effects are 
quite sizeable. In the peak region the bottom mass leads 
to a reduction and in the tail region to an enhancement 
resulting in a reduced negative slope of the distribution. 
As can be seen from Fig. [TTb . where the relative size 
of the mass effects with respect to the massless result 
is shown, the values of the thrust distributions are en- 
hanced up to 4% in the dijet tail region 0.05 < r < 0.15. 



to 0(a^) matrix elements and matching conditions and four-loop 
renormalization group evolution. 

For the soft function the dispersion method leads to an approx- 
imate result for the two-loop non-logarithmic corrections which 
is however sufficient for the brief analysis carried out here. 



In the peak region r < 0.05 the relative size of the re- 
duction of the distribution is even larger and leads to a 
sizeable effect. Although non-perturbative effects have 
to be included for a more definite numerical analysis, it 
is clear that the secondary bottom mass effects are im- 
portant for the analysis of event-shape distributions at 
lower cm. energies from JADE 30] or from current B- 
factory data. In the figure we have also indicated the 
field theory scenarios needed to describe the respective 
thrust ranges: for Q = 14 GeV and m& = 4.2 GeV only 
scenarios III and IV can arise because for very small r 
the bottom mass is located between the soft and the jet 
scale, which then increase for larger r. We have imple- 
mented the factorization theorem with a strict expansion 
of all two-loop corrections in the product of the matching 
conditions and the jet and soft functions, and there is no 
displacement of the blue solid curve at the transition from 
scenario III to scenario IV. The change of slope visible at 
the transition point of scenarios III and IV is related to 
the use of the different dynamical flavor numbers for the 
renormalization group evolution. 

An example where three scenarios are needed to de- 
scribe the full thrust spectrum is the case of secondary 
top quark production with mt = 175 GeV for Q = 



22 



500 GeV. Here the top mass is between the hard and 
the jet scale for small r, and also scenario II is required. 
The distributions with and without secondary top quark 
mass effects are displayed in Fig. [T2"r using the same la- 
beling as in Fig. [TT] We see that the behavior and the 
size of the secondary top quark mass corrections are very 
similar to the bottom quark case discussed before, albeit 
with the property that the peak region of the partonic 
distribution is confined to much smaller r values due to 
the large cm. energy Q. A more detailed numerical anal- 
ysis is presented in Ref. Q- 

VIII. CONCLUSIONS 

In this paper we have set up a method to deal with 
secondary massive quark radiation effects for inclusive 
jet cross sections. Depending on the mass value and the 
kinematic variables in the jet cross section the jet in- 
variant mass and the scale of soft large angle radiation 
can vary continuously crossing mass thresholds involv- 
ing the secondary quark pair. This makes the theoreti- 
cal method to deal with calculations of matrix elements, 
matching conditions and the summation of the correct 
large logarithmic terms complicated. We have set up 
an effective theory framework that combines the mass- 
less partonic modes described in soft-collinear effective 
theory with collinear and soft mass modes. While the 
collinear and soft modes for massless quarks and gluons 
in inclusive jets typically have different invariant masses, 
the collinear and soft mass modes have the same invariant 
masses. This leads to different effective theories depend- 
ing on the values of the mass and the kinematic variables. 
We have developed such an effective theory scheme for 
thrust in e + e~ collisions considering the effects of sec- 
ondary massive quark radiation (through gluon splitting) 



in massless quark pair production. The resulting scheme 
for thrust requires four different effective theories and al- 
lows for a continuous description of ultra-heavy quark 
masses in the decoupling limit down to the limit of van- 
ishing quark mass where the massless description is re- 
covered. The setup can thus be applied to every possible 
kinematic situation. 

The mass mode setup is in the spirit of the variable 
flavor number scheme originally developed by Aivazis, 
Collins, Olness and Tung where infrared-safe hard Wil- 
son coefficients were obtained with the subtraction of 
low-virtuality parton splitting terms to deal with quark 
masses in hadron collisions which are smaller than the 
hard scale but larger than Aqcd- Our method represents 
an effective field theory realization of the variable flavor 
number scheme using the framework of soft-collinear ef- 
fective theory for inclusive jet processes and its principles 
can be well applied to other processes involving heavy 
quark production including hadron collisions. 

For thrust and other observables where the invari- 
ant mass of the massive quark pair enters the observ- 
able, it is possible to describe the secondary massive 
quarks through a dispersion integral involving the ab- 
sorptive part of the massive quark vacuum polarization 
function and a massive "gluon" propagator. This allows 
to discuss the field theory setup for the simpler case of 
a spontaneously broken non-Abelian gauge theory com- 
bined with QCD involving massless quarks where the ad- 
ditional gauge bosons have a common mass. Conceptual 
as well as many technical aspects can be studied in great 
detail considering these massive "gluons" . In this work 
we have concentrated on developing the effective theory 
framework and the results concerning the massive gauge 
bosons, and we outlined the calculations required to ob- 
tain the results for secondary massive quarks. Other ap- 
plications shall be addressed in subsequent work. 
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